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Abstract 

We compute characteristic (quasinormal) frequencies corresponding to decay of 
a massive charged scalar field in a Reissner-Nordstrom black hole background. 

It proves that, contrary to the behavior at very late times, at the stage of 
quasinormal ringing the neutral perturbations will damp slower than the charged 
ones. In the limit of the extremal black hole the damping rate of charged and neutral 
perturbations coincides. Possible connection of this with the critical collapse in a 
massive scalar electrodynamics is discussed. 

Black hole (BH) quasinormal modes (QNM's) govern the decay of perturbations at 
intermediate times and are important when studying dynamics of black holes and external 
fields around them. Now, the interest in quasinormal modes is stipulated by three their 
features. The first is connected with possibility of observing of QNM's, and of obtaining 
the "footprint" of a black hole, with the help of a new generation of gravitational antennae 
which are under construction (see |l| for a review). The second concerns the Anti-de- 
Sitter/Conformal Field Theory (AdS/CFT) correspondence ||, from which it follows 
that a large black hole in AdS space corresponds to an approximately thermal states in 
the CFT. Thus the decay of the BH perturbations can be associated with the return 
to thermal equilibrium of the perturbed state in CFT 0. Note that in AdS space, 
quasinormal modes of black holes govern the decay of a field at late times as well, since 
there is no power-low tails and the decay is always exponential (see || and references 
therein). The third feature is related to possible connection of BH QNM's in some space- 
time geometries with the Choptuik scaling ||, 0, @, @. 

When Horowitz and Hubeny calculated massless scalar quasinormal modes for a four- 
dimensional SAdS black hole they found that for intermediate black holes, i.e. when the 
horizon radius is of order the Anti-de Sitter radius (see also 0), the striking approximate 
relation takes place 

w/ m r+, (1) 
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where r + is the horizon radius, 7 is the Choptuik scaling parameter for a black hole being 
formed from a collapse of a scalar field in free space. Then, exact connection between the 
Choptuik scaling parameter and the imaginary part of the QNF's was obtained for the 



massive scalar field in the background of the three-dimensional BTZ BH It gives us 
a hope that such a connection between characteristic frequencies of black holes, at least 
in special space-time geometries, and parameters describing situation near the threshold 
formation of a black hole exists for all black holes. Thus knowledge of such relations 
would be a most useful, since calculation of quasinormal modes is much simpler work 
then numerical simulation of gravitational collapse. 

Then it was observed in || that since the nearly extremal RN BH background is 
effectively described by AdS^ black holes under a spherically symmetrical dimensional 
reduction [H], the relation between the scalar QN modes of the nearly extremal RN BH 
and the Choptuik scaling may be obtained. 

Recently it has been found that at the stage of quasinormal ringing the neutral per- 
turbations, corresponding to a massless scalar field, will damp slower than the charged 
ones for RN, RNAdS, and dialton black holes ||. At very late times, on contrary, the 
charged perturbations will damp slower, thereby dominating in a signal as was shown in 
@. In addition the damping rates (imaginary part of u) of charged and neutral quasi- 
normal modes for massless scalar field coincide for the nearly extreme RN black hole 
0. This fact, although not understood, may be connected with the Choptuik scaling if 
remembering that the Choptuik scaling parameter 7 is the same (0.37) both for charged 
T(J and neutral [\ L| massless scalar field collapse. The case of a massive scalar field is 



more complicated from a critical collapse view. There both types of behavior appear (see 



12 1 and references therein), I (with a mass gap) and II (when the black hole being formed 
may be of infinitesimal mass) depending on ratio of the length scale of the initial data to 
the inverse Compton wavelength characteristic of the field. Yet the critical collapse was 
studied only for massless electrodynamics and for a massive scalar field. The complete 
picture of critical collapse for massive scalar electrodynamics is lacking as far as we are 
aware. This gives us one of the reasons for studying the quasinormal spectrum associated 
with decay of a charged massive scalar field. 

The other reason is to find out how a massive scalar field, interacting also electro- 
magnetically with a charged black hole, decay between an initial burst of radiation and 
the late time tails. In this direction the evolution of a neutral massive scalar field at 
late times was studied in |T3|] for a Schwarzshild background, in [[T7|], |Tj|, ]TS| for a RN 



background, and in jTJ] for the dilaton black hole. At the stage of quasinormal ringing, 



decay of neutral massive scalar field was studied in |2(J for Schwarzshild and Kerr black 
holes. The quasinormal modes of dilaton black holes were studied in 0, [[y|, |I6| . 

The complete picture of a charged massive scalar field evolution in a charged black 
hole background is lacking at either stages. 

We shall consider the evolution of the massive charged scalar field in the background 
of the Reissner-Nordstrom metric: 

ds 2 = -f(r)dt 2 + f-\r)dr 2 + r 2 d£ll (2) 

where f(r) — 1 — — + The wave equation of the complex scalar field has the form: 

<P, ab g ab - ieA a g ab (2 ( j>. b - xeA^) - ieA a , b g ab <P + /i 2 = 0, (3) 

here the electromagnetic potential A t = C — — , C is a constant. After representation of 
the charged scalar field into spherical harmonics the equation of motion takes the form: 

ip, tt + 2ie-i>, t - 4>, r * r * + Vijj = 0, (4) 
r 
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where 




the tortoise coordinate is defined by the relation dr* = jjrr, and tp = -0(r)e _,aJ ', uj = 
UR e — iuj m . The effective potential goes to a constant at the horizon and at infinity. We 
compute the quasinormal frequencies stipulated by the above potential using the third 
order WKB formula of S.Iyer and C.Will pTI: 



_ A( n ) - n(n) = n + - (6) 
2Q ' 2 



where A(n), Q(n) are second and third order WKB correction terms depending on the 
potential Q and its derivatives in the maximum. Here Q = —V + uj 2 — 2^-uj. Since Q 
generally depends on uj, the procedure of finding of the QN frequencies is the following: 
one fixes all the parameters of the QN frequency, namely, the multipole index I, the 
overtone number n, of the black hole Q and M, and of the field m and e; then one 
finds the value of r = r at which Q attains a maximum as a numerical function of uj and 
substituting it into the formula (J^j one finds uj which satisfies the equation (^J). Note that 
the effective potential, being frequency dependent, is, generally, complex. Nevertheless, 
following [[HJ , we treat uj as if it were real when finding ro and only continue Eq. (|6|) into 
the complex frequency plane. We deal here with n = modes as those dominating in a 
signal. 

First we examine the behaviour of the neutral massive scalar field in a RN back- 
ground. Recently in [|0| through numerical simulation of the wave equation solution it 
has been obtained that the relaxation process depends on the value M/i: when Mfj, <C 1 
the relaxation depends on the field parameters and does not depend on the spacetime 
parameters, while at Mfi ^> 1 the dependence on the black hole parameters appears. At 
the stage of quasinormal oscillations the evolution of the massive scalar field governed by 
QNF's depends both on black hole and field parameters. 

QNF's for different values of \i are presented on FigJT]- |3| for Z = 1,2,3. The I = modes 
are not given here, since the WKB accuracy is not sufficient in this case (nevertheless the 
main conclusions as to the quasinormal behaviour conserve in this case as well). The real 
part of the quasinormal frequency, i.e. the oscillation frequency, grows with increasing 
of the mass field fx, while the imaginary part of uj, representing the damping rate, falls 
down. Note that we are limited to the low-lying quasinormal frequencies whose real parts 
correspond to tunneling near the potential barrier (uj 2 fa V max ). In general, the maximum 
value of fj, (and of e for charged scalar field) depends upon the mode under consideration 

23. 

Next, we investigate the massive charged scalar field case. The real part of w as a 
function of imaginary part of uj is plotted for different / and Q for \i = 0.1 on Fig.4-6. 
We see that the neutral perturbations (e = 0) have less imaginary part than that of the 
charged perturbations. Thus the charged perturbations decay faster. This repeats the 
behaviour of the massless case 0. When increasing charge Q, uji m grows up to some 
maximum near Q = 0.8M and then falls down, while ujR e increases monotonically. The 
more e the less explicit this maximum. 

In addition for the near extremal black hole the imaginary parts of charged and neutral 
perturbations tend to coincide as in a massless case. In other words for the near extremal 
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Figure 1: RN QNF's for I = 1, n = 0, fi = 0,0.1,0.2,0.3, Q runs the values 
0,.1,.2,.3,.4, .5, .6, .7, .8, .9, .95, .99 




Figure 2: RN QNF's for I = 2, n = 0, fi = 0,0.1,0.2,0.3, Q runs the values 
0,.1,.2,.3,.4, .5, .6, .7, .8, .9, .95, .99 



black hole cj/ m does not depend on e for a charged massive scalar field. In this connection 
if one supposes some kind of relation between uji m of the nearly extremal RN BH and 
the parameter of the critical collapse in a massive scalar electrodynamics, then he could 
expect that ^massive being dependent on mass field n does not depend on e in the region 
where type II Jl2| critical behaviour occurs. Yet, at present we have much more questions 
then answers here. 
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Figure 4: RN QNF's for I = 1, n = 0, e = (bottom), 0.05, and e = 0.1 (top), Q runs 
the values 0, .1, .2, .3, .4, .5, .6, .7, .8, .9, .9, .95, .99, /x = 0.1 
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Figure 5: RN QNF's for I = 2, n = 0, e = (bottom), 0.05, and e = 0.1 (top), Q runs 
the values 0, .1, .2, .3, .4, .5, .6, .7, .8, .9, .95, .99, /i = 0.1 
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Figure 6: RN QNF's for I = 3, n = 0, e = (bottom), 0.05, and e = 0.1 (top), 
Q runs the values 0, .1, .2, .3, .4, .5, .6, .7, .8, .9, .95, 0.995, 0.999, ^ = 0.1. At Q = 0.999 
cu Im = 0.088517 at e = 0, uj Im = 0.088512 at e = 0.5, and u Im = 0.088470 at e = 0.1. 
Thus the coincidence of ujj m for different e in the near extremal regime is well within a 
third order WKB accuracy. 
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